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5jT)' Abstract. A general technique for the study of magnetic Rashba Hamiltonian in 

quantum graphs is presented. We use this technique to show how manipulating 
the magnetic and spin parameters can be used to create localized states in a certain 
periodic graph (T3 lattice). 

> 

g 1 Introduction 

In this work, we discuss the creation of eigenvalues in periodic quantum graphs by 
certain external interactions, namely, by magnetic field and spin-orbit coupling. 

The analysis of quantum graphs, i.e. of differential operators on singular 
one-dimensional manifolds, becomes one the central topics in the mathematical 
physics during last decades, see [28,29,31-33,35]. This has many reasons; in par- 
ticular, quantum graph Hamiltonians appear in the de Gennes-Alexander theory 
of superconductivity [4, 15,25]. Some other fields of applications are described 
e.g. in [31]. 

The spectral theory of compact quantum graphs has many common features 
with the usual theory of differential operators, cf. [6,42,48]. Nevertheless, such 
an analogy is rather limited when considering non-compact structures. Some par- 
ticular features of quantum graph models become obvious if one studies periodic 
configurations. For example, for a large class of periodic Schrodinger operators in 
Euclidian spaces the spectrum is known to be absolutely continuous [9,50], while 
even simplest periodic quantum graphs can have eigenvalues [12,33]. Some other 
examples may include the sensibility of periodic quantum graphs to some arith- 
metic characteristics [24]. 

Recently, in the physics literature one discussed the so-called extreme local- 
ization in the lattice (dice lattice) [1,54]. From the mathematical point of view, 
it was shown that under certain magnetic fields the Hamiltonian of a quantum 
graph with the corresponding shape has no bands of continuous spectrum, and 
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the spectrum consists of infinitely degenerate eigenvalues. This effect was ob- 
served also experimentally by transport measuring in superconducting and metal- 
lic wire networks [39-41]. (It is worth emphasizing that bound states in Tt, lattice 
appear without any external interactions [51]; the coexistence of the continuous 
and the point spectra is implied by the rich internal symmetry of the lattice and 
of its dual, Kagome [2].) Various aspects of this localization mechanism and 
its stability under disorder and external interactions are studied in a number of 
works [8,38,52-54]. In particular, it is shown that additional interactions, like the 
inter-particle interaction, destroy the extreme localization mechanism, and con- 
tinuous spectrum appears [53]. 

In Ref. [7] it was shown that in some periodic quantum graphs similar lo- 
calization phenomena can be induced not only by magnetic fields, but also by 
spin-orbit interaction at certain values of the Rashba constant [14,47]. Never- 
theless, this analogy is limited, and the numerical analysis of [8] shows that the 
Rashba localization does not appear in the Tt, lattice. In the present paper we 
consider the above situation with both the spin-orbit and magnetic interactions. 
We note that the quantum graph models with spin were studied previously e.g. 
in [10, 11, 13,27], but the attention was mostly concentrated on Dirac- and Pauli- 
type operators. The theory of Rashba Hamiltonians is not developed even in the 
Euclidian spaces, where the spin-orbit interaction promises to show effects which 
are absent in the scalar case, like embedded eigenvalues in short-range potentials 
or localization in crystals, giving hence possibilities for constructing new nanode- 
vices [17,34]. 

The aim of the paper is two-fold. First, we are going to describe the 
Schrodinger operators in two-dimensional networks with magnetic field and spin- 
orbit interaction. An essential part here is the reduction of the quantum graph 
Hamiltonian to a certain discrete equation. In the scalar case, an analogous proce- 
dure was done in [23] for the solutions of the stationary Schrodinger equation and 
recently in [44] for the spectra. Note that there is another approach to the relation- 
ship between the quantum graphs and tight-binding Hamiltonians coming from 
some asymptotic considerations [36,37]. Second, by considering the localization 
effects in the -lattice we would like to attract the attention of researchers work- 
ing on quantum graphs to potential applications in the study of superconducting 
networks. In section|2]we give a mathematical formalism of quantum graphs with 
external interactions; essentially we describe rigorously the constructions of the 
works [7, 8]. In section[3]we use this machinery to study the spectrum of the T$- 
lattice with a magnetic field and the Rashba interaction. We show that the spectral 
problem is of supersymmetric type and that the study of some energy levels is 
equivalent to the study of zero modes in a certain discrete model. As a result, we 
give a rigorous justification of the extreme localization for the case of non-trivial 
scalar potentials on the edges and non-ideal couplings at the nodes. We show that 
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at zero spin-orbit interaction this effect is independent of the edge potential. At the 
same time, it appears that the generic Rashba interaction destroys the localization. 
We also show that at a certain combination of the magnetic and spin parameters a 
new localization effect appears, where one can localize one of the spin projections 
using the magnetic field. 

2 Quantum graphs with external interactions 

2.1 Schrodinger operator on a quantum graph embedded in 
Euclidian space 

In this section we describe the construction of the Hamiltonian in a two- 
dimensional network with magnetic field and Rashba interaction. Recall that the 
Rashba Hamiltonian of a two-dimensional system acts on two-component vector 
functions as takes the form [14,47] 

H= (p-A) 2 + 2^(p-A,crxn)+[/, (1) 

where A is the magnetic vector potential, kn is the Rashba constant expressing 
the strength of the spin-orbit interaction, U is a scalar potential, cr is the vector of 
Pauli matrices, and n is the unit vector orthogonal to the plane of the system. The 
second term, which is the formal mixed product, on the right-hand side of CQ) takes 
into account the spin-orbit coupling. For Ur — the problem splits in two identical 
scalar problems. The corresponding Hamiltonian for a network is obtained by 
projecting all the interactions onto each edge and by introducing suitable boundary 
conditions at the nodes, which will be described below. (We remark that some 
effects of the Rashba interaction and the magnetic field in a wire can be studied in 
other types of models [19].) 

Let "V be a uniformly discrete subset of the xy-plane in R 3 , the set of nodes 
(vertices). The uniform discreteness means the existence of a constant d > such 
that | a - (3\ > d for all a, (3 G "V with ol^ (3. Denote 

lu/3 '■= \ol-(3\, e a f3 ■= j—(f3-oi). 

l af3 

Some nodes are connected by a directed edge. The set of all edges will be denoted 
by S, S C "V x "V . The edge with initial vertex a. G "V and terminal vertex f3 G "V 
will be denoted by a/3. For at E ~f denote indega := #{(3a G $ }, outdegc* = 
#{ol(3 G $}, deg ex := indeg o: + outdeg ex. We assume that the degrees satisfy the 
following conditions: that 

there exists N G Z with 1 < deg a < N for all aeY; (2) 
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in particular, we assume that there are no isolated vertices. The configuration 
consisting of all segments [a, (3], a/3 G £ will be referred to as a metric graph or 
as a (wire) network. We assume that the system has no self-intersections and that 

< mf{l a p} < sup{/ Q/3 } < oo. (3) 

The quantum state space corresponding to the metric graph is introduced as 
follows. Each edge a/3 will be identified with the segment [0,/ a «] such that a 
is identified with and (3 is identified with l a p. The state space of each edge 
a/3 is Jrf? a p := L 2 ([0,Z a/ g],C 2 ). The state space of the whole structure is Jt? = 

On each edge consider a real- valued scalar potential U a p G L 2 [0,l a p]. To 
avoid unnecessary technical difficulties we will assume that the scalar potentials 
are uniformly L 2 -bounded, 

sup||£/ a/ 3|| L 2 < °o. (4) 

Assume that the system is subjected to an external magnetic field given by a vec- 
tor potential A G C^M 3 ,]^ 3 ). This induces magnetic potentials on each edge, 
a a p(t) ■= (A(a + te a p),e a p). 

Denote by k R the Rashba constant. The spin-orbit interaction can be taken into 

account by adding the term 2k R (^ — i— —a a p(t)j (cr x n,e aj g) with n = (0,0, 1). 

Therefore, the dynamics along each edge a/3 is described by the differential ex- 
pression 



L a p = +2k R (-i—-a a p(t)^((Txn,e al 3)+U a p 

= (i-£+a a p{t)+k R o a p) +U a p-k 2 R , 

where 

( e a p 2 + ie a p X \ 

p \e a 02-ie a pi J 

For a uniform magnetic field with the strength B G R 3 it is useful to use the sym- 
metric gauge, A(r) = - B x r. In this case the magnetic potentials a a p are con- 
stant, a a p = ^(Bxa,e a p). 

Denote by L an operator in Jt? acting as 

(fa/3) >-»■ (L a pf a p) (5) 
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on functions i a p G H 2 ([0, l a p], C 2 ) satisfying at each qgT: 



fafl(0) 



f 7 a(Z 7a ) =:f(a), a/3, 7aG#, 



E 

apes 



dt 



-i(a a p + k R o a p) i a p(0) 



Y, -^-'A a i3a + k R op a ) fp a (lp a ) = e(a)f(a), 



PaeS 



(6a) 
(6b) 



where e(a) are real- valued parameters. The case £(a) = may be considered as 
an ideal coupling, which is the analogue of the Kirchhoff coupling in the scalar 
case. We are going to consider L as the Hamiltonian of the system, and our next 
aim is to show its self-adjointness. 



2.2 Self-adjointness and spectral analysis 

Denote by Qi the set of all functions f = (f a /3), with components (f a p) £ 
Q)H 2 ([0,l a p],C 2 ), a/3 G <§, which are continuous at all nodes, i.e. such that 
the condition (l6al) is satisfied. Clearly, for f G ^ the values f(a), a 6/, have the 
direct sense. Furthermore, for f G & and ct G Y denote 



f'(Q) 



E 

apes' 



Idt 



i(a a p+k R o a p) f a p(0) 



- 52 ~u ~K a /3a +k R Op a ) fp a (lp c 



Pa.eS 



Consider in a linear operator fl with domain acting by the rule ©. 

Proposition 1. The operator fl is closed. For any f G domll = 3) the vectors Ti : = 
(f(a)) and T'f : = (f (a)) belong to £ 2 (Y,C 2 ), and the map (T,P) : domL - 
f(Y , C 2 ) © £ 2 (f , C 2 ) w surjective. For any f, g G domn fftere /ioM? 



(f,n g ) - (nf, g ) = (rf,r g ) - (rf,r g ). 

Proof. Denote by aj g the unitary transformation of J^ a p given by 
©a/3f(0 = exp (i jf (a a/3 (» +*RCT ai g)ds)f(0. 



(7) 



(8) 



Denoting d := — we see (5 -ia a/9 -ik R o a p)® a p = & a pd. 
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By the Sobolev inequality, for any c\ > there exists C2 > such that for any 
/ > and <p G # 2 [0, /] there holds 

||(p'|| 0O <C 1 / 1 /2|| (p -|| L2[a/] + _^|| (?) || L2[a/] . 

Note that for any * G [0,/ a/3 ] one has ||f a /3(0llc 2 = ll®a/3 f a/3(0 lie 2 - Therefore, 
using the above estimate, for any i a p G // 2 ( [0, l a p] , C 2 ) one has 

||fa/3(0||c2 = ||®a/3 f a/3(')|lc 2 

^ c l'a| 1 1 d 2 ©a/3 f «/3 I U Q/3 + TT^ 1 1 ®a/3 f «/3 1 1 

'a/3 

= C l / a / Jll a/3 - ia <*/3 ~ '^ R O a p) 2 i a p \\ ^ + || a/3 fa/3 1 1 ^ 

'a/3 

= C 1 Z a/3 1 1 ( ^ ~ ifl a/3 - ^a/3 ) 2 fa/3 1 1 + 1 1 fa/3 1 1 ,^ q/3 

'a/3 

and, in the same way, 

\\(d —M a p — ikR(7 a p)f a p(t)\\ C 2 

^ C ^\\i d -Map-ikROapftapWje^ + 4^ || f a/3 ||^ • 

'a/3 

Using the assumptions © and © we conclude that there exist positive constants 
C\ and C2 such that for any a/3 G S, f a p G // 2 ( [0, l a p] , C 2 ) , t G [0, Z a/ g] one has 

\\l a p{t)\\<C l \\L a pf a p\\+C 2 \\f a pl (9a) 
||(a-ia a/3 -i^(j a/3 )f Q/3 (f)|| <Ci||L aj3 / aj g||+C2||/ ai a||. (9b) 

Here the norms are taken in C 2 on the left-hand side and in ^ a p on the right-hand 
side. 

Denote by II the operator acting in by the rule © on the domain domll = 
©a/3e<f-^ 2 ([0^a/3];C 2 ). Clearly, II is closed. By (l9aT). the linear maps 

T a p 7 : domfi 3 f ^f a ^(0) -f 7 «(/ 7 a) e C 2 , a/3, 7a G 

are bounded with respect to the graph norm of II. Therefore, the restriction of II 
to the subspace where all these functionals vanish is a closed operator. As this 
restriction is exactly II, the operator II is closed. 
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For f £ 3 the inclusions IT, T'f £ £ 2 (^, C 2 ) follow immediately from the 
estimates © and the assumption ©, and the identity © can be verified directly 
using the partial integration. 

To prove the surjectivity condition, we fix first four functions /)& £ # 2 [0, 1] 

with/Jty) = 5^4/, i,y,Jfe,/ £ {0, 1}. Take arbitrary £ £ 2 (^,C 2 ). Denote 



r a/3 := exp (ijf " (fl aj g(j) +^C7 a/3 )d^ 

= exp(ij ^ a a p(s)ds^cosk R l a/3 +io a/3 smk R l a( ^) e\J(2). (10) 

By direct calculation, the function f £ J^ 3 whose components are of the form 
fa/3 = ©a/3ga/3, where 



+ ^/io(7^)^ , («)-^/ii( 7 ^)^ , (/3) ) 
dega \l a/3 J deg/3 \Z a/3 / ap 



'—/ml 

f a/3 / aegp \f a/3 - 

lies in and satisfies (IT, T'f) = (£, £') . □ 

Proposition [l] shows that the space 8f := £ 2 (^,C 2 ) and the maps r,r' : 
domll — > 5f form a boundary triple for II, see e.g. [21,46] for a detailed discus- 
sion. The self-adjointness of II would follow from the following assertion [21]: if 
II has at least one self-adjoint restriction (i.e. if II* is symmetric) and A is a self- 
adjoint operator in Sf , then the restriction of II to the vectors (p £ domll satisfying 
T'<p= AT<p is self-adjoint in Jff. 

Consider the restriction D of II to the functions f satisfying IT = 0. Clearly, 
this restriction is nothing but the direct sum a ^ £( f D a p, where D a p is an opera- 
tor in Jf a p acting as f a/3 h-> L af3 { a/3 on functions satisfying f a/3 (0) = l a p{l a p) = 
0. As each D a p is self-adjoint, so is D. Note that L itself is the restriction of II to 
the functions f satisfying T'f = 7Tf, T = diag (e(cx)) . This implies 

Proposition 2. The spin-orbit Hamiltonian L is self-adjoint. 

To carry out the spectral analysis of L it is useful to relate the resolvents of L 
and D by Krein's resolvent formula [21], 

(D-E)- 1 - (L-E)- 1 = r(E) [M{E)-T\- X y*(E) (11) 

where E £ specLU specD and the maps y(E) and M(E) axe defined as follows. 
For a given E specD and £ £ 1 2 {Y , C 2 ), the function y(E)£ = (f a/ a) is the solu- 
tion to (n -£)f = satisfying IT=£. The map M(E) : £ 2 (r,C 2 ) -> £ 2 (y,C 2 ) 
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is given by M(E) = T'y(E). A direct consequence of Eq. (fTTT) is the relationship 

specL\spec£> = ^E £ specD : G spec (M(E) - T) j. (12) 

Moreover, E specD is an eigenvalue of L iff is an eigenvalue of M(E) — T, 
and y(E) ker (M(E) — T) is the corresponding eigensubspace [26]. 

Denote by D a p the self-adjoint operator in L 2 [0, l a p] acting as g a p i— ► — g^Vg + 
~ kfygap on functions g ai g G // 2 [0,/ aj g] satisfying the Dirichlet boundary 
condition, g a p(0) = g a p(l a p)- Note that the operators ©Jj,#D a /3®a/3 are of the 
form D a p ® D a p- I n particular, the spectra of D a p coincide with those of D a p 

and are discrete sets, and specD = U a /3e<f s P ec ^a/3- 

Eq. (fT2l) shows that the spectrum of L outside specD is completely described 
in terms of M(E). The question whether specD or some parts of it enter to the 
spectrum of L must be analyzed individually taking into account the magnetic and 
spin parameters and the topological properties of the graph. 

Therefore, to carry out the spectral analysis for L it is useful to calculate the 
map M(E). This can be done in terms of special (scalar) solutions to the equation 

-y" + U a py = zy, zeC. (13) 

Namely, denote by s a p and c a p the uniquely determined solutions of (fT3l satis- 
fying the boundary conditions 

s a p (0; z) = c' a/3 (0; z) = 0, s' af3 (0; z) = c a p (0; z) = 1 . 

Now let £ G £ 2 (y,C 2 )- To find y(£)£ =: (f a/3 ) we need to solve the boundary 
value problems 



(i^ + a ap + kRO af3 ) 2 + U a p - k\ 



f a p=Ef a p, 



f«j9(0) = e(«), f a p(l a i3)=&P)- 

Writing f aj g := © a /3g a /3> where & a p is the unitary transformation from (EJ), we 
rewrite (TUT) as a boundary value problem for g a p, 

-g'Lp + u <xpg a p = (E + ki)g a p, 
g a p(0) = |(a), g a p(l a p) = v^((3) 

The solution to <TT~5l) takes the form 
* a/3 (?;£ + k\) 



(15) 



s a p(l a p;E + kl) 



? a p€(P)-c a p(l a p\E + k 2 R )£(a) 

+ c afj {t-E + k z R )i{a). 
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Now we have 
g«/3(0) = 

g' a p( l ap) 
Noting that 

we arrive at 



1 

s a p(lap',E + k R ) 

1 

s ap(l<xp',E +k R ) 



T&ptiP) ~ Ccxp(Lp;E + 4)£(a) 

J a p(i<*p;E+k 2 R )T* af3 £(p)-Z(cx 



(16) 
(17) 



f(«)= £ g«/3- £ -cp a g' a p{i a p) 

aP<ES PaeS" 



<pm 



xP^S s ap(laP',E + k R 



+ £ n p I 7 2 \ T /3a^(/3) 

PaegSpaVPat'Z + K R) 

c a p(l a p;E + kl) ^ s'^lp^E + kl) 



+ 1 



£(«)• (18) 



*/3e<f s ap{lap\E + k R ) p ae g S/3 a (lp a ;E + k R )- 
Using Krein's resolvent formula (fTTl) we come to 

Theorem 3. The set specL\ specD consists exactly of the real numbers E such 
that e spec [M(E) - T], where M(E) and T are operators in f{Y, C 2 ), M(E) 
is given by (fT51) and T = diag (e(oc)). Moreover, such E is an eigenvalue ofL iff 
is an eigenvalue ofM(E) — T, and y(E) ker (M(E) — T) is the corresponding 
eigenspace. 

We remark that in the above calculations it does not matter whether £ is in 
£ 2 or not. Actually, all the construction hold for any set of vectors £(«) G C 2 , 
ex. E y . This observation can be formulated as follows: 

Theorem 4. For E ^ specD, any continuous solution f to (L — E)l = has the 
form 



MO 



s a p{t-,E + kl) r m (/ £ + £2) f (a) " 

s a p(l a p;E + k R ) L ^ J 

+ 0a/3C a/3 (?;£ + 4)f(«)- 



Swc/z a function satisfies the boundary conditions (l6bl z/f 

aPeg^OLpV-ap^^-Kjt) C 



tfle^ s a.p\ i a.p^ + K R) Po.eS s Pa\ i Pa^ + K Rl 



v c a(3 (l af3 -E + kl) v s'^jlp^E + kl) i 

£ O ,r I l2\ L 77 p , 7 2 N + £ ( Q: ) f («)- 

t p e g s ap\ l aP> t! ' -tKr) f3 ae g SpuKlpaiC' -r Kr) j 
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Such an f is an eigenfunction ofL (i.e. belongs to L 2 ) iff (f(a)) G £ 2 (y, C 2 ). 

Note that similar formulas for more simple situations were obtained earlier 
e.g.in [3,5,8,22,23,32]. 

The expression (fT8l) can be simplified significantly if all the edges are the 
same, i.e. if l a p = I and U a p = U, s a p = s, c a p = c for all ol(3 G S. Note that in 
this case the spectrum of D coincides with the Dirichlet spectrum of the operator 
D := —A 2 /At 2 + U — k\ on the segment [0, /] and hence is a discrete set. We have 



M{E)£(a) 



s(l;E + k 2 ) 

outdega c (I ; E + kl) +mdeg as' '(/; E + k\) £(ct)\. (19) 



Even this expression admits further simplifications. 

Proposition 5. Assume that all edges are identical, l a p = I, U a p = U,U is even, 
U(t) =U (l—x), and the coupling constants £(a) are of the form £(a) = deg a E, 
then spec L\ spec D = t E l (spec A), where t £ (E) = c(l;E + k 2 i ) + £s(l;E + k 2 ,) and 
A is the discrete Hamiltonian, 



acting on the space t Li (f ', C 2 ;deg) with the scalar product 



(£,»7>deg= £ dega-£(a) 77(a). 

Proof. If the potential U is even, one has c(l;E-\-kfy =s'(l;E + kf,) :=t(E), see 
e.g. [44], hence 



M(E)-T 



1 



A - * e (£) deg , deg = diag ( deg a) , 



s(l;E + k 2 ) L 
where A is the discrete Hamiltonian in ir(Y, C 2 ), 
&(«) = £ £ 

The condition G spec \M(E) — T 7 ] takes the form G spec [A — t E (E) deg] in 
£ 2 (r,C 2 ), which is equivalent to 0g spec [A -t E (E)] in £ 2 (r,C 2 ; deg). □ 
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Figure 1: A finite piece of T3 -lattice 



Proposition 13 shows that the spectral problem for a class of quantum graphs 
reduces to the study of the tight-binding Hamiltonian A. In the case U = 

and £(<*) = one has t £ (E) = cos E + kj^, and we arrive at specL = 

Arccos 2 spec A — k\ (up to the discrete set specD), which is exactly the formula 
connecting the network and the tight-binding spectra in the de Gennes-Alexander 
model of superconductivity [4]. For scalar situation, an analogue of this corre- 
spondence was given e.g. in [5] for the Laplacian on compact graphs, in [16] for 
the Laplacian on non-compact graphs, and in [44] for more general Schrodinger 
operators. At the same time, proposition |5] does not exhaust all possibilities of 
such a reduction, i.e. the reduction to a discrete Hamiltonian is possible also for 
some non-even U. (Such questions will be discussed in greater detail in [45].) 
One of such situations will be discussed in the next section. 



3 Spectrum of T3 -lattice 
3.1 Description of the lattice 

In this section, we consider the spectral problem for a quantum graph whose un- 
derlying structure is the so-called ^-lattice (see figured)- The nodes are the points 

.3 v3 . .3 v3 . 

Ot m , n , Pm,n, lm,n With CL m ,n = mSL\ + na 2 , *\ = (-> Y ^ * 2 = ^2 1 ~2~ 
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1 V3 

/3 m ,n = a m)B + (1,0,0), 7 m ,„ = a m; „ + (-, — ,0), m,n e Z, i.e. 

_ /3(m + n) v^n-ra) \ /3(ra + n) + 2 v^n-ra) \ 

am '""l 2 ' 2 ' /' P m > n -{ 2 ' 2 ' /' 

'3(m + «) + l y/3(n — m+l) 



"fm.n 



-,o). 



The edges are 



€m,n,4 — G-m,nfitn,n—\i e m,n,5 = < ^-m 1 n^lm+\,n—\i e m,n,6 ~ CX-m,nfim,n- 

All the edges have the length 1. The direction vectors of e mM j are e 7 = 
(cos y,siny,0), j = 1,...,6. 

3.2 Reduction to tight-binding Hamiltonian 

We will assume that the system is subjected to the following external interactions. 
On each edge there is the same potential U G L 2 [0, 1]. The lattice is subjected 
to the uniform magnetic field B = (0, 0, 2%E, ) orthogonal to the plane, and the 
magnetic vector potential in the symmetric gauge is A(x) = (—n%X2, K^xi,0). In 



what follows we use the magnetic parameter co 



expressing the magnetic 



flux through the elementary rhombus (for example, CKm,«/3m,nO:m,n+ i7mrc)- 

The external magnetic field implies non-trivial magnetic potentials on e m ^ n j, 

a mA i = <o(2ra + n), a m ^ 2 = Co(m + 2n), a m ^ 3 = co(n-m), 
a m ,n,4 = ~(0{2m + n), a m ^ 5 = -a(m + 2n), a m ^ 3 = Co(m- n). 

The dynamics along e mM j is described by the differential expression 

>,2 



d 2 

L m ,n,j = 0-~j^ a m,n,j "H O m ,n,j) ~\~ U fc^, 



m,n,j 



ej2 + ieji 
e j2 -ieji 



exp 







exp 



3 2- 
m,nGZ, j e {1,...,6}. 



i(--^) 
. K 2 3 ' 
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Here are the components of the vectors ey and k R is the Rashba constant. 
We consider boundary conditions of the form © at all the nodes assuming that 
there are only two types of coupling constants: X := e(oc m>n ) and /i := e(/3 ra ,„) = 

The corresponding matrices T OT)W j from Eq. (TTOb . T w ,nj = exp [i(a w , n j + 
kRO m ,n,j)] , are as follows: 

T m , M = <, nj4 = e ira ( 2 ^ ^ }, (20a) 

^,2 = C,5 = e lffl(m+2 ' ,) {cos^ + isin^( e ^ /6 6 ^ }, (20b) 

^,3 = C )6 = e^-'"){cos^ + isin^( el ^ /2 6 ^ }. (20c) 

Clearly, for any m, n G Z one has 

outdega w>n = 6, indeg/9 m ,„ = indeg7 m ,„ = 3, 
indega,„,„ = outdeg/3 m ,„ = outdeg7 OT) „ = 0. 

For the subsequent analysis we use the fact that the lattice is bipartite. Rep- 
resent the set of nodes as the disjoint union "f = "Pq U Y\, "Pq = {ac m n }, "V\ = 
{flm,n} U {lm.n}- Clearly, for the set of edges one has S <zYqxY\. With re- 
spect the the decomposition f(Y : C 2 ) = £ 2 (%,C 2 ) ®P(V\X 2 ) the operator T 
in theorem|3]takes the block-diagonal form, 

X 
/I 

Using the above decomposition and Eq. ( fT9l) . we rewrite M(E) — T as 

a(E) = 6c(l;E + k 2 R ) + Xs(l;E + k 2 R ), 
b(E) = 3s\l;E + k 2 R )+Lis(UE + k 2 R ). 



where 



A5.\~ s (rn,n) = tm,n,lf(cX-m,n) tm,n+l, 3^(^/7?, n+1 ) ~\~ ^m— \,n-\-\^X\^-m— l,n-(-l)) 

A*f(a m , n ) =T; ; „ ;1 f(7 m ,n) + < )n)2 f(/3 m _i )B ) + < )n)3 f( 7m , n _i) (21) 
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The operator A*, A* : £ 2 (y u C 2 ) -> £ 2 {%X 2 ) is adjoint to A, A : £ 2 (%X 2 ) -> 
£ 2 ('^ / i,C 2 ). Using theorem |3] we write the condition £ e specL or, equivalently, 
OG spec (M(E)-T), as 



Note that E specD in all the above constructions, where D is the Dirichlet re- 
alization of — d 2 /dt 2 + U — kjt on [0, 1]. The question whether specD is a part of 
specL or not admits a simple answer in our case. 

Lemma 6. For all (0 and one has spec D C spec L. 

Proof. Using the Schnol-type arguments, cf. [12, 33], it is sufficient to show that 
for each E E specD the equation Li = Ei has a bounded solution f satisfying the 
boundary conditions ©. Let E E D and g be the corresponding eigenfunction of 
D. Choose any infinite path without intersection on the graph and any non-zero 
vector z e C 2 . For a fixed e E £P set f e := ® e gz, where & e is given by (JU). Now 
extend f to the whole graph in such a way that (a) f e = for e & and (b) on each 
b E 8? one has f b = ® b gz b , where the vectors z b are chosen in such a way that the 
boundary conditions © are satisfied. By construction, there holds Li = Ei. At 
the same time, due to the unitarity of the matrices T m)W j the obtained function f is 
bounded. This finishes the proof. □ 

3.3 Supersymmetric analysis 

Eq. d22l is a typical supersymmetric spectral problem. Using propositionfTTI and 
corollary[l2]in Appendix one easily sees that the set E of E for which the condition 
(|22|) is satisfied is the union E = Ei U E2 U £3, 




(22) 



E <£ specD : a(E)b(E) ^ and a(E)b(E) E specA*A, 




To summarize, 



Proposition 7. specL = Ei U E2 U E3 U specD. 
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Note that the sets E2, £3, and specD are discrete. Therefore, only the set X4 is 
responsible for the continuous spectrum. Writing f(a m; „) =: f(m,n), we note that 
A* A is an operator on £ 2 (Z 2 , C 2 ) of the form 

A*Ai(m,n) =6f(m,n) 

+ (*m,n,l*m,n+\3 + ^nfi^n+l,^ (m,n + 1) 

+ (^m,n,3^m,n-l,l + ^m,n,4^m,n-l,6)f(m,n- 1) 

+ ( T m,n,5 T m+l,n,3 + T m,n,6 Wl,n,2) f ( m + 1,*) 

+ ( T m,«,2 T m-l,n,6 + ^,„,3 T m-l,n,5) f (*» ~ 1 , «) 

+ ( T m,n, 1 tm- 1 ,n+ 1 ,5 + ?m,n,2 X m- 1 ,n+ 1 ,4 ) f (m - 1 , W + 1 ) 

+ ( T m,n,4 W 1 ,n- 1 ,2 + )B)5 T m+ 1 ;M _ i , i ) f (m + 1 , n - 1 ) , 



i.e. 



A*A = 6 + cos ft) • sh^fc/j • A 
+ 2 cos ft)- cos 2 ]cr — sin 2 



COS (ft) -|) 







cos(ft) + -) 



A 
A 



where A is a spinless operator in £ 2 (Z 2 ), 

A/(m, n) = e- 3i(0m f(m, n+l) + e 3i(0m f(m, n-l) 
+ e 3l0)n f(m + l,n)+ ^ n f{m - 1 , n) 
+ e - 3ifi) ( m+n )/(m - 1 , n + 1) + e 3iffl ( m +")/(m + 1 , n - 1 ) , 

and 

Af (m, n) =7?i e - 3ifflm f (m, n + 1 ) + R\ e 3iam f (m, n - 1 ) 
+7? 2 e 3iffl "f (m + 1 , n) + R* 2 e- 3i03n f(m - 1 , n) 
+R 3 e - 3iffl ( m +") f (m - 1 , „ + 1 ) + R* e3i»(m+«) f( m+ 1) 



with 



#1 



V 



o 

3 y/3. 



3 
2 





#3 = 



#2 



/ 



-V3i 



( 

V 

V3i 






3 >/3. 
2 + 



3 , v^.\ 
2 + ^ 







/ 



(23) 



The expression for A* A shows explicitly the contribution of the magnetic and spin- 
orbit parameters to the spectrum. Let us discuss the situations where the spectrum 
shows certain localization phenomena. 
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Figure 2: Classification of the nodes (3 m ,„ and 7 m , n . 

3.4 Magnetic field induced extreme localization 

If the spin-orbit interaction is not taken into account, Icr = 0, then one has 

K 

A* A = 6 + 2 cos (O (A © A) . In particular, at (0 - — G %Z one has A* A = 6, i.e. 

71 

the spectrum of A*A degenerates to a point. If CD — — G TtZ but the spin-orbit 

interaction is non-trivial, similar phenomena occur only at certain values of the 
Rashba constant, i.e. Icr G TtZ. For generic values of kR, obviously, there are some 
bands of continuous spectrum. 

71 

Let us analyze the sets £2 an d £3 for this case, i.e. for CO — — £ TtZ and 

kR G ttZ. Clearly, the set E2 is empty, as ^ specAM. Let us look at the operator 
AA*. 

71 

Lemma 8. For co — — e TtZ and kR e TtZ one has e specAA*. 

Proof. In view of periodicity, it is sufficient to show that the equation AA* y = 
has non-trivial bounded solutions, ty/ G £°°(^i;C 2 ). 

Let us classify the nodes f3 m ,n and j mM as shown in figure |2j Consider all 
vector-valued functions on Y\ vanishing at the white marked nodes. For such a 
function i/a, the condition A*\y — is of very simple form, because in the expres- 
sion <f2Tb only two of the six terms on the right-hand side are non-zero. Therefore, 

fixing the value of i// at a single black marked node one uniquely extends \y to a 

71 

bounded solution of A* i// = 0. The conditions co G — + %7L and kR G kZ guarantee 

that this solution is well defined, i.e. that the phase factor along each cycle on the 
hexagonal lattice of black nodes is 1 . □ 
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To summarize the previous considerations we note that the set specD consists 
of the real E satisfying s(l;E + k R ) = 0. Proposition[7]reads as 
% 

Theorem 9. Let CO £ — + TiZ and k R £ itL, then the spectrum ofL consists of the 
real numbers E satisfying at least one of the following conditions: 

fc(UE + kl) + ^s(UE + e R )y^XUE + 4) + ^s{UE + kl)^ =i (24a) 

s\l;E + k 2 ) + ^s(UE + k 2 )=0, (24b) 

s(l;E + k 2 R )=0. (24c) 
Each point of the spectrum is an infinitely degenerate eigenvalue. 

Note that the eigenvalues (I24al) are the most interesting ones, as they arise as 
the limit of the continuous spectrum. The Dirichlet eigenvalues d24cl) are usually 
neglected in the physical works. 

In the simplest case, when the scalar potential is zero and the couplings are 

trivial, i.e. U = 0, X = fl = 0, one has s(x;E) = —= sin y/Ex, c(x;E) = cos y/Ex, 

. . VE 

and Eq. (I24al) takes the form 

COS 2 yjE+k 2 R E {0 ,1,1}, 

which was previously obtained in [54] for the case k R = 0. 

We remark that the presence of the extreme localization is periodic with re- 
spect to the shifts k R \-^ k R + n, but not the energy levels themselves, as the func- 
tions s(-;E + k R ) etc. are not periodic with respect to the Rashba constant. It is 
worthwhile to note that the above results hold for any potential U and any coupling 
constants A and /i. 

3.5 Magneto-spin induced localization 

Another interesting situation appears at cosk R = 0, i.e. at k R £ — + nZ. In this 
case one has 

cos (co-—) A 

3 71 

cos(co + -)A^ 

For the values CO £ — — + TtZ the first component of A* A degenerates. In 

6 

particular, any function of the form (/,0), / £ £ 2 (1?) becomes an eigenfunction 

of A* A. For co £ — + 7tZ the same holds for the functions (0,/). 
6 

For further analysis we calculate the spectrum of A. 
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K 

Lemma 10. For CO G ±— + 7rZ ?/ze spectrum of A is absolutely continuous and 

6 

covers the segments [ — 2 a/3, —a/3] and [y/3, 2 a/3] . 
Proof. Consider the unitary transformation 

U : f(Z 2 ) 3 (f(m,n)) ^ (e 3ia>nm f(m,n)) G £ 2 (Z 2 ). 
Clearly, U is unitary, and the operator A := U*AU has the form 
A/(m, n) =e- 6imm f(m, n + 1 ) + e 6imm f(m, n - 1 ) + f(m + 1 , n) + f(m - 1 , n) 



+ e 



-3i«(2m-l 



/(m-l,n+l)+e 



3i«(2m+l 



/(m+l,n-l). 



The operator obtained has the same spectrum as A, but is periodic with respect 
to the shifts n 1— > n + 1 and can be studied using the Bloch analysis. Making the 
Bloch substitution f{m, n) — e lw? ^ m , where q G [0, 2k) is the quasimomentum, we 
observe that the spectrum of A is the union of the spectra of operators H(q) acting 
in £ 2 (Z) and defined by 

(#(«)£) m = 2exp[-3ifi>(in-i)+i|] cos [3o(m-i)-|] £ m _i 
+ 2cos [600m — q\ E, m + 2exp 3ift)(m + ^) — i^ cos 3tt)(m + ^) — ^ £m+i- 

The operators are nothing but the Harper operators for the triangular lat- 

K 

tice [18]. Note that for CO G ±— + ftZ all these operators are invariant under the 

6 

shift m\-^m + 2. Therefore, substituting into the equation H{q)E, = Et, a vector t, 
satisfying t, m = e 10 <^ m _2 for all m, where G [0,2k) is another quasimomentum, 
one arrives at a 2 x 2 linear system for the components £0 and £1, 



,i(3ffl/2+?/2) 



— + -)e ^i+cos^o 



+ e 



W2-*/2)cos(^-|)& 



e-W-^) CO s(^-|)^o-cos^ 1 



+ e -i(3«/2+,/2) cos (^ + 1)^=^ 



The condition for the determinants to vanish takes the form 

,2 



£ 2 



— = l+cos q -\- cos(^ — 0) cos ^ 

3 1 2 

= - + (cosg + -cos(#- 0)) + sin(q-0). 

Taking here all possible values of q and we arrive at the conclusion. 



□ 
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LemmaflOlmeans that for values 0) in question, the spectrum of A* A has a con- 
tinuous part, which is the union of the segments [0, 3] and [9, 12], and an infinitely 
degenerate eigenvalue 6. 

Therefore, we arrive, as in subsection 13 .41 to a series of infinitely degenerate 
eigenvalues E satisfying the same equation (I24al) (i.e. the same eigenvalues as in 
the extreme localization case), which are isolated in the spectrum, but we have 
additionally bands of continuous spectrum given by 

(^l^ + ^ + ^^l^ + ^j-^l^ + ^ + l^l^ + ^GtO^JUt^]. 

In particular, for the free case with zero coupling constants one has the following 
characterization for E to be in the spectrum of L: 

cos 2 sJe^ r e [0, 1 -] u { \ } u [I |] U { 1 } . 

The localization effect described in this subsection seems to be not covered by 
the existing works, and it would be interesting to know whether it can be really 
observed. As for different values of the magnetic parameteres we have completely 
different eigensubspaces of A* A, we conjecture that this localization mechanism 
can be used to control the spin polarization by the magnetic field, but this needs a 
further analysis. 



4 Appendix. Supersymmetric spectral analysis 

Here we prove the following proposition. 

Proposition 11. Let M% are some Hilbert spaces, A be a bounded linear 
operator from M{ to Mj,, and m £ R. On 3tf{ © ffii consider the operator 

m A 
^A —m 

Then 



specL = — J spec(AA* +m 2 ) U y spec(AM + m 2 ). (25) 

This proposition is formulated (without proof) in [43] and is nothing but an 
abstract version of proposition 2.5 in [49]; we give here a complete proof just for 
the sake of completeness. 

Proof. First note that specAA* \ {0} = specAM \ {0} [20]. Clearly, 



Ll =(""«"" 2)- (26) 



A*A + m 2 
AA* + m' 
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Therefore, specL 2 \ {±m} = spec(AA* +m 2 ) \ {m 2 }, and for any A G specAA* \ 
{0} = specAA* \ {0} at least one of the numbers — VA +m 2 , \fX-\-m 2 lies in 
specL. Let us show that actually they both are in the spectrum of L. 

Let A > 0, A G spec A*A, then there exist a sequence (0„), § n G M'x such that 
||0„|| > 1 and lim(A*A - X)§„ = 0. Denote 



Wn 







(27) 



Clearly, 



A* 
A 



which implies 

By direct calculation, 



I Vwll > ^||0n|| > 



(28) 



(L - v^) Wn = ( a/AW - m) ( (AM " A } M 



Therefore, lim(L — VA + m 2 ) = 0. Together with d28t this implies VA +m 2 G 
specL. 



To show — vA+m 2 G specL one has to consider the functions 



Vn := 



A-( V / A + 



m 2 — mj 



A* 
A 



where ||0 n || > 1 and lim(AA* — A)0„ = and to repeat the above steps. To finish 
the proof of Eq. (l25l) it is necessary to study the points ±ra. 

For m = 0, Eq. (l2"6l) reads as specL 2 = specAA* U spec A* A, and the conditions 
G specL and G specAA* U spec A* A are equivalent. 

Assume m ^ and m G specL, then there exist sequences (0 n ) G J^i, (<p M ) G 

with 

||<M + ||<Pn||>l (29) 



and 



lim(L — m) 



0ii 



lim 



A*<p ;j 



0. 



(30) 



jp n J \A<j) n - 2m(p n 

Clearly, this implies limA*A0,, = 0. Assume that lim^)„ = 0, then (l30l) shows 
lim<p„ = which contradicts (l29l) . Therefore, there exists a subsequence (<j>' n ) of 
(0 n ) such that || > e for some e > 0. Together with limA*A0/ = this implies 
G specA*A. 
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Assume now G specA*A, then there is a sequence (0 n ) G ^4 with ||0 n || > 1 
and lim(A*A0„, n ) = lim ||A0„ || = 0. Then 

lim(L-„,)(*")=li m ( A J=0, 

from which m G specL. 

The relationship between the conditions — m G L and G spec AA* can be 
proved in a completely similar way. □ 

It may be useful to have an alternative formulation of proposition [TT] 

Corollary 12. There holds 



specL\ {— m,m} = — y spec(AA* + m 2 ) U y spec(AA* +m 2 ) \ {— m,m} 
= —\J spec(A*A + m 2 ) U \/ spec(A*A + m 2 ) \ {— m,m}. 

Furthermore, for m^O one /zas; m G specL (JO 6 specAM, — m G specL iff 
G spec AA*, and for m = ?/zere /zoMs G specL zJO G spec A* A U spec AA*. 
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